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Abstract 

We investigate the average characteristic polynomial E [ YiiLi ~ )] where 
the real random variables which form a determinantal point process as- 

sociated to a bounded projection operator. For a subclass of point processes, 
which contains Orthogonal Polynomial Ensembles and Multiple Orthogonal 
Polynomial Ensembles, we provide a sufficient condition for its limiting zero 
distribution to match with the limiting distribution of the random variables, 
almost surely, as N goes to infinity. Moreover, such a condition turns out to 
be sufficient to strengthen the mean convergence to the almost sure one for 
the moments of the empirical measure associated to the determinantal point 
process, a fact of independent interest. As an application, we obtain from a 
theorem of Kuijlaars and Van Assche a unified way to describe the almost sure 
convergence for classical Orthogonal Polynomial Ensembles. As another appli- 
cation, we obtain from Voiculescu's theorems the limiting zero distribution for 
multiple Hermite and multiple Laguerre polynomials, expressed in terms of free 
convolutions of classical distributions with atomic measures. 



1 Introduction and statement of the results 
1.1 Introduction 

For any A^> 1, letxi,...,XArbea collection of real random variables which forms a 
determinantal point process associated with a rank N bounded projection operator. 
This means there exists for each a Borel measure /x^v on M and a /i^r (g) /x 7v-square 
integrable function Ki\f : M x M — M such that the joint probability distribution on 
of xi, . . . ,xn reads 

1 - ^ 



KN{xi,Xj) Y\^^N{dx.i), (1.1) 

together with the fact that the operator acting on L'^{fiN) by 

TTN ■■ f ^ J KNix,y)f{y)nN{dy) (1-2) 
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is a (non- necessarily orthogonal) projection on an A^-dimensional subspace of L^(/i7v)• 
Strictly speaking, this is not a standard way to introduce determinantal point 
processes, but it is easy to obtain from the standard references on the subject 
[1, 29, 31, 45] that a determinantal point process (in the usual sense) with a kernel 
satisfying the latter conditions induces such random variables, and vice versa. To 
these random variables, we associate their average characteristic polynomial. 



Xn{z) = E 



N 

Z Xj 

i=l 



zeC, (1.3) 



where the expectation E refers to (1.1), and we ask the following question: What 
is a sufficient condition so that the asymptotic distribution of the zeros of xn and 
the limiting distribution of the random variables Xj's coincide as — t- oo ? More 
precisely, if one denotes by zi, . . . , z^v the (non-necessarily real nor distinct) zeros of 
Xn and introduces the zero counting probability measure 



N 

i 

UN 



1 ^ 



N 

i=l 



the purpose of this work is to investigate the relation between the weak conver- 
gence of vn and the almost sure weak convergence of the empirical measure of the 
determinantal point process, namely 



Let us first observe that if xi, . . . , x^v are i.i.d real random variables, say, with law 
fiN and (finite) mean niN, then this relation is trivial. Indeed, note that u^- = Smj^ 
since the Vieta's formulas yield xn{z) = {z — ttt-tv)^- Thus un and almost surely 
fi^ converge towards the same limiting probability distribution as — >• oo if and 
only if TTi-Tv converges to some real number m and fi^ almost surely converges weakly 
towards 5m- If /^7V has no atoms, then the Xj's form a determinantal point process, 
with kernel Kfyf{x,y) = 6xy (the Kronecker's delta) and reference measure /^at, that 
is associated to the trivial projection on the 0-dimensional subspace of L?'{iin)- We 
emphasize this trivial case does not satisfy our hypotheses and will not be further 
discussed. 

The situation is much more interesting in the case of non-trivial determinantal 
point processes. It is for example known that the eigenvalues of an iV x random 
matrix drawn from the Gaussian Unitary Ensemble (GUE) is a determinantal point 
process, and that xn is the A^-th monic (i.e. with leading coefficient one) Her- 
mite polynomial. After an appropriate rescaling, the zero distribution vj^j converges 
weakly towards the semi-circle distribution as A^ — )■ oo, and so is almost surely the 
spectral measure fi^ , although the numerous proofs of these two facts seem quite 
independent. 

In this work we provide a sufficient condition so that, as N ^ oo, the convergence 
of the moments of z^at is equivalent to the almost sure convergence of the moments 
of fi^ for a large class of determinantal point processes, see Theorem 1.3. We will 
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actually show that condition implies the simultaneous moment convergence of z^at 
and of the mean distribution E[/i^], defined by E[/i'^](A) = E[/i'^(j4)] for any 
Borel set ^ C M, and moreover forces the moments of jl^ to concentrate around 
their means at a rate A^^"*"*", see Theorem 1.8. At this level of generality, the latter 
concentration result is new and may be of independent interest. To do so, we 
develop a moment method for determinantal point processes, involving weighted 
lattice paths, see Section 2.2. 

Besides the theoretical aspect, and as we shall illustrate in Section 3 and 4, 
such a statement provides two useful practical consequences. On the one hand, 
the almost sure convergence investigation for such determinantal point processes is 
thus reducible to the asymptotic analysis for the zeros of polynomials, for which 
analytic tools have been developed in special cases. On the other hand, one can use 
the probabilistic background of the random models to obtain a description of the 
limiting zero distribution of average characteristic polynomials, which in particular 
cases happen to be special functions of interest in other areas of mathematics. 

Before stating our results, let us first introduce and discuss what is already known 
for two important classes of point processes that will be covered by this work. 

1.2 Orthogonal Polynomial Ensembles 

Examples of Orthogonal Polynomial (OP) Ensembles are provided by eigenvalue 
distributions of unitary invariant Hermitian random matrices, including the GUE, 
Wishart and Jacobi matrix models; they also arise from non-intersecting diffusion 
processes starting and ending at the origin. In the latter examples, has a density 
with respect to the Lebesgue measure. They moreover play a key role in the res- 
olution of several problems from asymptotic combinatorics, such as the problem of 
the longest increasing subsequence of a random permutation, the shape distribution 
of large Young diagrams, the random tilings of an Aztec diamond (resp. hexagone) 
with dominos (resp. rhombuses). This time /i at is a discrete measure. For further 
information, see [35, 32, 33] and references therein. 

The joint distribution of real random variables xi, . . . ,XAr drawn from an OP 
Ensemble can be written as 

1 2 ^ 

^ ^<i<i<N 1=1 

where Zj^ is a positive normalization constant and ^at is a measure on M having all 
its moments. One can rewrite that distribution in the form (1.1) by introducing the 
symmetric kernel 

N-l 

KN{x,y) = ^ Pk,N{x)Pk,N{y), (1-6) 
fc=o 

where Pk,N is the /c-th orthonormal polynomial for /iTv- This is the kernel associated 
with the orthogonal projection onto the subspace of L^(/xjv) of polynomials having 
degree at most — 1. 

An important observation, provided by a classical integral representation for OPs 
attributed to Heine, see e.g. [16, Proposition 3.8], is that the average characteristic 
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polynomial XN associated to an OP Ensemble equals the A^-th monic OP with 
respect to /xat- Since OPs are known to have real zeros, vn is thus supported on M. 

As we shall recall in Section 2, the mean distribution E[/i^] of a determinan- 
tal point process equals jj'KN{x,x)fj,]\[{dx). Quite remarkably, it turns out that 
when Kn has the form (1.6), the convergence of the mean distribution has been 
investigated in the approximation theory literature, where it is referred as the weak 
convergence of the Christoffel-Darboux kernel. Using the determinantal point pro- 
cesses terminology, Nevai [42] and Van Assche [46] actually proved that the simul- 
taneous weak convergence of E[/i^] and holds for OP Ensembles as soon as a 
growth condition on the recurrence coefficients of the -Pfc^iv's is satisfied (a defini- 
tion for the recurrence coefficients of OPs is provided in Section 3). Nevertheless, 
their proofs involve the Gaussian quadrature associated to OPs, an argument which 
does not seem to be generalizable to more general determinantal point processes. 
More recently, and in the case where the supports of the measures fi^ are uniformly 
bounded, Simon [43] proved the simultaneous moment convergence of E [/i^] and 
by means of elegant operator-theoretic arguments, which have been of inspiration 
for this work. 

Concerning the almost sure convergence, after ordering the Xj's and Zj's, Dette 
and Imhof [19] were able to obtain in the case of the CUE, that is the OP Ensemble 
associated to fi]y{dx) = e~^^^/^dx, an upper bound for P(max^-|^ \ xi — Zi\ > e), from 
which the almost sure simultaneous convergence of fi^ and z/jy follows. They estab- 
lished that the same picture holds for the Wishart matrices, that is for iij\[{dx) = 
x^"e~^^l[o^+oo)(2;)dx where a > 0. Dette and Nagel [20] also worked out the Jacobi 
case, namely /iAr(dx) = {1 — x)^°'{l + x)^^l[_i^i]{x)dx, where a,f3 > —1. Their re- 
sults moreover cover the associated /3-Ensembles. In both works, the proofs strongly 
use the explicit tridiagonal matrix representation for these random matrix models, 
which is unfortunately not available for more general unitary invariant Hermitian 
random matrix models nor determinantal point processes. 

Let us now describe a larger class of determinantal point process that will be 
also covered by this work. 



1.3 Multiple Orthogonal Polynomial Ensembles 

Firstly introduced by Bleher and Kuijlaars [9] to describe the eigenvalue distribution 
of an additive perturbation of the GUE, breaking the unitary invariance, Multiple 
Orthogonal Polynomial (MOP) Ensembles show up in several perturbed matrix mod- 
els [8, 10, 18], in multi-matrix models [38, 22, 23] as well, and in non-intersecting 
diffusion processes with arbitrary prescribed starting points and ending at the origin 
[39]. For general presentations, see [36, 37] and the references therein. 

The joint distribution of real random variables xi, . . . , xat distributed according 
to a MOP Ensemble has the following form 



(xj — Xij det 



i<i<i<Af 



\x) Wi^N{xj)\ 
^ J i,j=l 



\ x) Wr,N{Xj) > 

J i,j=l 



N 



HfiNidxi), (1.7) 



i=l 
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where //at is a measure on M having all its moments, Zn^^r is a normalization constant 
and the weights wi^n, ■ ■ ■ ,Wr,N G L'^{iJ'n) are such that (1.7) is indeed a probability 
distribution. The multi-index n = {ni, . . . ,nr) G depends on and satisfies 
l^i=i iT'i = where throughout this paper we denote N = {0, 1, 2, . . .}. Note that 
we recover OP Ensembles by taking r = 1. 

It turns out, see Section 4.5, that there exists a sequence (i-*fc,Af)fceN of monic 
polynomials with degP^^Ar = k, and a sequence {Qk,N)keN of (non-necessarily poly- 
nomial) L^(/ijv )-functions which are biorthogonal, that is 

j Pk,Nix)Qm,Nix)fiN{dx) = Skm, /c, m € N, (1.8) 

such that we can rewrite (1.7) in the form (1.1) with the kernel 

N-l 

KN{x,y) = ^ Pk,N{x)Qk,N{y)- (1.9) 

k=0 

Kuijlaars [36, Proposition 2.2] established that the average characteristic poly- 
nomial xn associated to (1.7) is the n-th (type II) MOP associated with the weights 
Wi,N, 1 < i < and the measure ^at, see Definition 4.1. 

The simultaneous convergence of the empirical measure fi^ and the zero distri- 
bution i/AT of the associated MOPs is expected for several MOP Ensembles. It is 
for example the case for non-intersecting squared Bessel paths with positive starting 
point and ending at the origin. Indeed, for this MOP Ensemble E[/i^] converges 
towards a limiting measure described in terms of the solution of a vector equilib- 
rium problem, see [39, Theorem 2.4 and Appendix], and the limit of Vf^ benefits 
from the same description [40]. The same situation holds in the two-matrix model 
with quartic/quadratic potentials, by combining the works [22] and [21]. For the 
non-intersecting squared Bessel paths model, which is equivalent to a non-centered 
complex Wishart matrix model, the almost sure convergence of jl^ towards the so- 
lution of the vector equilibrium problem has moreover been obtained from a large 
deviation principle in [28]. For the two matrix model, to prove a large deviation 
upper bound involving a rate function associated to a vector equilibrium problem is 
still an open problem, see [24] for further discussion. For these two determinantal 
point processes, the almost sure simultaneous convergence of fi^ and will be a 
consequence of what follows, see Remark 1.9. 

It is now time to describe the general setting for which our results hold. 

1.4 Statement of the results 

Consider a collection of real random variables xi, . . . ,xn which forms a determi- 
nantal point process associated with a rank N bounded projection operator ttat on 
L?'{iin), with kernel ITat. Under these only assumptions, Ktv is defined ® ^n- 
almost everywhere, which is not sufficient to characterize a determinantal point 
process, i.e. (1.1) is not well defined. We then follow [29, Remark 5] and generalize 
their approach in the following way. The spectral decomposition for compact opera- 
tors [44, Theorem 1.4] provides two biorthogonal families {Pk,N)k=o and {Qk,N)^=o 
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of L'^{fj,j\f), namely which satisfy 

Pk,N{x)Qm,N{x)flN{dx) = 6km, < k,m < N - 1, (1.10) 
such that the following equality holds in L'^{fiN /^Af) 

KN{x,y) = ^ Pk,N{x)Qk,N{y)- (1-11) 



k=0 



We take the right hand-side of (1.11) as our definition for Kj\f. Note that although 
the -Pfc.Af's and Qk,N^s are still only defined /iAr-almost everywhere, the probability 
distribution (1.1) now reads 



1^ r 1^ A 

Pfc_i,Ar(xi) det Qk-i,N{xi) m(dxi) (1.12) 

J i,k=l I J i,k=l 

1=1 



and is properly defined. Thus, our class of determinantal point processes matches 
with the Biorthogonal Ensembles introduced by Borodin [11], where we emphasize 
that the Pfc.Af's and the Qk,N^s are L^(;U7v)-functions. 

Now, given a sequence of such determinantal point processes indexed by N (the 
number of particles), 

Pk,N , Qk,N , fJ-N? , (1-13) 

J k=0 V J k=0 

K ) N>1 

we assume moreover the following structural assumption to hold. 
Assumption 1.1. 

(a) For each N , the two families {Pk,N)k=o {Qk,N)k=o can be completed in 
two infinite biorthogonal families {Pk.N)^^^ and {Qk,N)k^f^ °f L'^ifJ'N), that is 
which satisfy 

Pk,N{x)Qm,Nix)lJ-N{dx) = 6km, A;, 771 G N. (1.14) 

(b) There exists a sequence (qAr)Ar>i of integers having sub-power growth, that is 
for every n > 1, 

qjv = o(A^^/") as iV ^ oo, (1.15) 

such that for all A; G N, 

xPk,N{x) G Span Pm,N{x) 

V / m=0 

Remark 1.2. OP and MOP Ensembles both satisfy Assumption 1.1 (with = 1 
for all > 1). A class of determinantal point processes which satisfy this assumption 
but for which q^v may grow is provided by mixed-type MOP Ensembles, originally 
introduced by Daems and Kuijlaars to describe non-intersecting Brownian bridges 
with arbitrary starting and ending points [15]. Delvaux showed that the average 
characteristic polynomial XN is in this case a mixture of MOPs [17]. 
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Let P be the probability measure associated to the product probabihty space 
(g)^(M^,Pjv), where (M^,Piv) is the probabihty space induced by (1.12). The 
central theorem of this work is the following. 

Theorem 1.3. Assume there exists e > such that for every n > 1, 



max 

fc.mGN; |4-l|<e, l?f-l|<e 



as N ^ oo. Then, for all £ &N, 



xPk,N{x)Qm,N{x)fiN{dx) 



o(7Vi/") (1.16) 



lim 



/.'A"(d.,-/.V(d.) 



0, ¥- almost surely. (l-l^) 



In practice, the sub-power growth condition (1.16) may be interpreted as the 
condition that a strong enough normalization for the Xj's has been performed. 

Remark 1.4. Assumption 1.1 (a) and (b) provide together for each N the unique 
decomposition 

fc+qiv 

xPk,N{x) = {xPk,N,Qm,N)L^(^,j.f)Pm,N{x), A) G N. (1.18) 

m=0 

Thus (1.16) is a growth condition for the coefficients lying in a specific window of 
the infinite matrix (i.e. operator on ^^(N)) associated to the operator f{x) i— )• xf{x) 
acting on Span(Pfc,Af)A:eN- 

As announced in the introduction, let us now provide more precise statements 
concerning the concrete uses of Theorem 1.3. Having in mind that probability mea- 
sures on M with compact support are characterized by their moments, the following 
consequence of Theorem 1.3 may be of use to obtain almost sure convergence results. 

Corollary 1.5. Under the assumption of Theorem 1.3, if there exists a probability 
measure fi* on M characterized by its moments such that for all i G 'N, 

lim / x^i']\f{dx) = / x'''fj*{dx), 
N^ooJ J 

then ¥ -almost surely fi^ converges weakly towards fi* as N ^ cc. 

As an example of application, we will obtain from a result of Kuijlaars and 
Van Assche a unified way to describe the almost sure convergence of classical OP 
Ensembles in Section 3, see Theorem 3.1. 

Similarly, when one is interested in the limiting zero distribution of xn, the 
following corollary will be of help. 

Corollary 1.6. Under the assumption of Theorem 1.3, if 
(a) for all N large enough xn has real zeros. 
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(b) there exists a probability measure fi* on M characterized by its moments such 
that for all ^ G N, 



then vn converges weakly towards fi* as N ^ oo. 

As an example of application, we will obtain in Section 4 a description for the 
limiting zero distribution of multiple Hermite and multiple Laguerre polynomials, 
see Theorems 4.5 and 4.6. At the best knowledge of the author, this is the first time 
that a description of these zero limiting distributions is provided in such a level of 
generality. 

Remark 1.7. Although it is not hard to see from our proofs that Theorem 1.3 
continues to hold for determinantal point processes on C (with the introduction of 
complex conjugations where needed), Corollaries 1.5 and 1.6 are not true in the 
complex setting. Indeed, consider the eigenvalues of an x unitary matrix dis- 
tributed according to the Haar measure, which are known to form an OP Ensemble 
on the unit circle with respect to its uniform measure. We have xn{z) = , and 
thus vn = 6o for all A^, but the spectral measure jl^ is known to converge towards 
the uniform distribution on the unit circle as A^ — )• oo. 

On the road to establish Theorem 1.3, we prove the following inequality which 
basically allows to extend the mean convergence of the moments of fi^ to the almost 
sure one. 

Theorem 1.8. Under the assumptions of Theorem 1.3, for every < a < 1 and 
any £ G N, there exists Ca/ independent of N such that for all 5 > 0, 



If moreover qjv and the left-hand side of (1.16) are bounded (seen as sequences of 
the parameter N), then (1.20) also holds for a = 1. 

Remark 1.9. Having in mind real OP and MOP Ensembles, let us stress that 
our results combine nicely with a Deift-Zhou steepest descent analysis. Indeed, it 
is known for such ensembles that one can represent in terms of the solution 
of a Riemann-Hilbert problem, see [16] (resp. [36]) for OP (resp. MOP) Ensem- 
bles. This, in principle, allows to use the Deift-Zhou steepest descent method, which 
yields a precise asymptotic description of ATat, and related quantities. In particular, 
the (icPfc^TVi Qm,iv)L2(/ijv)'^' '^hich turn out to be the recurrence coefficients of OPs 
and MOPs, see Sections 3 and 4, can be expressed in terms of the solution of the 
Riemann-Hilbert problem (see [16, (3.31)], resp. [26, Section 5]) and a control of 
their growth would follow from that steepest descent analysis. Such an asymptotic 
analysis also typically provides the locally uniform convergence and tail estimates 
for Ki^ as N ^ oo, from which would follow (1.19), and where the limiting measure 
^* has in general compact support. In most cases, the zeros of xn are real; this 




(1.19) 




(1.20) 
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is always true for OPs and also for important subclasses of MOPs, like Angelesco 
or AT systems. Thus, if one assumes the latter to be true, the combination of a 
successful Deift-Zhou steepest descent analysis together with Corollary 1.6 and The- 
orem 1.8 would provide the almost sure weak convergence of the empirical measure 
fi^ , and moreover the weak convergence of the zero distribution of the (M)OPs 
towards n* , without extra effort. For example, the determinantal point processes 
associated to the non-intersecting squared Bessel paths and the two-matrix model 
with quartic/quadratic potential discussed in Section 1.3 both satisfy the latter; the 
asymptotics of the recurrence coefficients are actually explicitly described in [40, 
(1.11)] and [21, Theorem 5.2] respectively. An other example of MOP Ensemble 
where the recurrence coefficients are explicit is provided by [5], in relation with the 
six- vertex model. 

The rest of this work is structured as follows. In Section 2, we establish Theorems 
1.3 and 1.8. In Section 3, by combining Corollary 1.5 and a result concerning 
the zero convergence of OPs obtained by Kuijlaars and Van Assche, we provide a 
unified description for the almost sure convergence of classical OP Ensembles. In 
Section 4, after a quick introduction to MOPs, we use Corollary 1.6 and Voiculescu's 
theorems in order to identify the limiting zero distribution of the multiple Hermite 
and multiple Laguerre polynomials in terms of free convolutions, and moreover derive 
algebraic equations for their Cauchy-Stieltjes transform. 



2 Proof of the main theorems 

In a first step to establish Theorems 1.3 and 1.8, we express all the quantities of 
interest in terms of traces of appropriate operators. 



2.1 Step 1 : Tracial representations 

Consider a determinantal point process associated to the rank N bounded projector 
TTjv acting on L'^{^n) with kernel Kat given by (1.11), so that 



N-l 



Im(7rAr) = Span Pk^N , Ker(7rAr) = Span Qk,N , 

\ / fc=0 \ / fc=0 

The usual definition of a determinantal point process, see e.g. [31], provides for any 
n > 1 and any Borel function / : M" — )• M the identity 



E 



r 1 " TT 

/(xi, . . . ,x„) det KN{xi,Xj) LuAr(dxj), (2.1 

1=1 



where the summation concerns all pairwise distinct indices taken from {1, . . . , A^}. 
Let M be the operator acting on L?'{^n) by 



Mf{x) = xf{x). 



(2.2) 



Then the following holds. 



9 



Lemma 2.1. For any i gN, 



E 



N 



Tr(7rArM^7rAr). 



Proof. By using (2.1) with n = 1, (1-11) and the biorthogonahty relations (1.10), 
we obtain 



E 



■ TV 

.1=1 



N-l „ 

/ X^Pk,N{x)Qk,N{x)fJ.N{dx) 



k=0 

N-l 



'TrN)Pk,N, Qk,N)L^ 



{^J■N) 



k=0 



Tr (vTArM^Ar). 



□ 



We also represent the variance of the moments in a similar fashion. 
Lemma 2.2. For any ^ G N, 

1 



Yar 

Proof. We write 
Yar 



x^fi^idx) 



N 
.1=1 



E 



Ar2 



■ N 

.1=1 



Tr(7r^M2Vjv) - Tr(7r7vMV7vM%iv 



,2£ 



+ E 



E 



Xj^Xj 



E 



■ AT 

E^ 

,i=l 



in order to obtain, thanks to (2.1) with n = 2 and Lemma 2.1, 

N 



Yar 



E 



J=i 

Finally, observe that 



Tr(7rivM^^7rAr) - JJ x^y^KNix,y)KNiy,x)nNidx)fiN{dy). 



X y KN{x,y)KNiy,x)fiN{dx)nNidy) 

N-l 



^ / ( / ^^'^^^y^y^^^^^^y^^^'^'^y^ Qk,N{x)^N{dx) 



k=0 
N-l 



k=0 

to complete the proof. 



□ 



We now check that the average characteristic polynomial xn equals the charac- 
teristic polynomial of the operator tt^Mttn acting on Im(7r7v)- 
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Proposition 2.3. //det stands for the determinant of endomorphisms oflm{Trj\f), 
then 

Xn{z) = det (z — ttnMttn), z G C. 



Proof. On the one hand, Vieta's formulas provide 



E 



N 

4 = 1 



Z Xi 



N 



n=l 



E 



il 7^ ■■■ 7^i„ 



and (2.1) provides for any 1 < n < 



E 



E 

«1 ¥= ■■■¥=^n 



det 



1 " TT 

K{xi,Xj) fiN{dxi). 

J «,.7 = 1 . 



1=1 



On the other hand, since tt^Mttn is an integral operator acting on Im(7rjv) with 
kernel {x,y) i— )• yKiy{x,y), the Fredholm's expansion, see e.g. [27], reads 



N 



det {z - ttnMttn) = + Y1 "i 



N-n 



n=l 



det 



1 " TT 

XjKN{Xi,Xj) /XAr(dXj), 

-I « J=l -T 



i=l 



from which Proposition 2.3 follows. □ 

The next immediate corollary will be of important use in what follows. 
Corollary 2.4. For any £ G N, 

xVjv(dx) = ^Tr((7r,vM^jv)0- 
The second step is to rewrite the traces in terms of weighted lattice paths. 
2.2 Step 2 : Lattice paths representations 

We introduce for each the oriented graph Qn = {Vn,£n) having Vn = for 
vertices and for edges 

£n = ^{n,k) ^ {n+l,m), n, /c G N, < m < /c + qTvj. 

To each edge is associated a weight 

WN(^{n,k) (n+l,m)^ = {xPk,N,Qm,N) L^i^j.^^), 

and the weight of a finite length oriented path 7 on is defined as the product of 
the weights of the edges contained in 7, namely 



(2.3) 



Then the following holds. 
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Lemma 2.5. For any £ G N, 

E x'fi^dx) =-Yl E ^^(^)' (2.4) 

fc=0 7:(0,A:)-^{^,fc) 

where the rightmost summation concerns all the oriented paths on starting from 
{0,k) and ending at {£,k). 

Proof. It follows inductively on i from (1.18) and the definition (2.3) that 

{TTNM^7TN)Pk,N ( ^^(t) ) ^'m.TV, ^, € N. (2.5) 

m=0 \7:(0,fc)-5>(£,m) / 

Thus, we obtain from the biorthogonality relations (1.10) 

TY(7r7vMV7v) = J2{{7rNM^TTN)Pk,N,Qk,N)L^ii,j^) 
k=0 

N-1 

= E E ^^(^)' (2-6) 

fc=0 7:(0,A:)-5-(£,fc) 

and Lemma 2.5 follows from Lemma 2.1. 

□ 

Next, we introduce 

Dn = |(n,m) E : m > ivj (2.7) 

and obtain a similar representation for the moments of v^. 
Lemma 2.6. For any £ G N, 

hr'uM{dx) = -Yl E ^^W- (2.8) 

fc=0 7:(0,fc)-s>(^,A:), 7ni:>jv=0 

Proof. Similarly than for (2.5), we have 

N-lf \ 
(tTnM ■ ■ ■ TTNM^7rN)Pk,N = Y] | Yl WNh)\Pm,N, /c G N. 



m=0 \'y:{0,k)^{e,m),'yr\Dff^ 



Since tt'L = ttn, this yields 



(2.9) 



N-1 

Tr((^jvM7r^)^) = ^ {{ ttnM ■ ■ ■ ■KNM nN)Pk,N, Qk,N)LH^^^) 

k=0 ^ } ' 

N-1 

= E E ^^(^) (2-10) 

fc=0 7:(0,fc)-5-(^,fc), 7nD]v=0 

and thus Lemma 2.6, because of Corollary 2.4. □ 
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If we denote by 7(m) the ordinate of a path 7 at abscissa m, then we can 
represent the variance of the moments of fi^ in a similar fashion. 



Lemma 2.7. For any i £ 



Yar 



1 iV-l 

k=0 7:{0,fc)-5-(2£,fc), 'y{i)>N 



E 



WAT (7). 



(2.11) 



Proof. We have already shown in (2.6) that 

N-l 



fc=0 7;(0,fc)-!>(2^,fc) 



(2.12) 



Since 



m=0 \ 7;(0,fc)^(f,m), ■y(t)<N 



we moreover obtain 



Tr(7rivM^7rArMVAr) = ^ ^ ?yiv(7). 

A:=0 7:(0,A:)-5-{2^,fc), ■y(l)<N 

Lemma 2.7 is then a consequence of Lemma 2.2 and (2.12)-(2.13). 

We are now in position to complete the proofs of Theorems 1.3 and 1.8. 

2.3 Step 3 : Majorations and conclusions 

Let us first provide a proof for Theorem 1.3 assuming that Theorem 1.8 holds. 
Proof of Theorem 1.3. It is enough to prove that for any given ^ G N 



(2.13) 



□ 



lim 



E 



X vn{Ax] 



0, 



(2.14) 



since (1.17) would then follow from Theorem 1.8 and the Borel-Cantelli lemma. As 
a consequence of Lemmas 2.5 and 2.6, we obtain 



E 



N-l 



I x^/i^(dx) - / xV^(dx) = 1 ^ ^^(7). (2.15) 



Since by following an edge of ^at one increases the ordinate by at most the 
rightmost sum of (2.15) will bring null contribution if k is strictly less that A^ — 
Observe moreover that the vertices explored by any path 7 going from (0, k) to (£, k) 
for some N — < k < N — 1 such that 7 n D^v 7^ form a subset of 

^{n,m) : 0<n<e, N - < m < N + CiN^y 
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As a consequence, if one roughly bounds from above the number of such paths by 
{2q]\f£Y, one obtains from (2.15) that 



E 



x'''i']\f{dx) 



< max 



I \^ 

\{xPk,N,Qm,N)L^{f,r^)\ ■ (2-16) 



It then follows from (2.16) together with the growth assumptions (1.15) and (1.16) 
that (2.14) holds, and the proof of Theorem 1.3 is therefore complete up to the proof 
of Theorem 1.8. □ 

We now prove Theorem 1.8 by using similar arguments than in the proof of 
Theorem 1.3. 

Proof of Theorem 1.8. Again, because following an edge of Qm increases the ordinate 
of at most qAT, the rightmost sum of (2.11) brings zero contribution except when 
k > N — c|7v£. Observe also that the vertices explored by any path 7 going from 
(0, k) to {2£, k) for some — < k < N ~ I and satisfying > N form a 
subset of 

|(n,m)GN2: 0<n<2i, N - Iq^i < m< N + 2qNiy 
As a consequence, we obtain from Lemma 2.7 the (rough) upper-bound 



Yar 



< 



x^fi^idx) 



max 



|(xPfc,7V,Qm,Af)L2(/.iv)f^- (2.17) 



Using the sub-power growth/boundedness assumptions on q^ and on the left-hand 
side of (1.16), Theorem 1.8 then follows from (2.17) and the Chebyshev inequality. 

□ 



3 Applications to OP Ensembles 

Consider a sequence of real OP Ensembles, introduced in Section 1.2, namely a 
sequence of determinantal point processes of the type (1.13) where the measures /iat 
on M have infinite support, have all their moments, and where Pk,N = Qk,N stands 
for the A;-th orthonormal polynomial with respect to ^n. We recall that it satisfies 
Assumption 1.1 with q^v = 1 and that xn and Pn.n have the same zeros. For every 
A^, the celebrated three term recurrence relation for orthonormal polynomials reads 

xPk,Nix) = ak+i,NPk+i,Nix) + bk^NPk,Nix) + ak,NPk-i,N{x), k>l, 
xPo,Nix) = ai^NPi,Nix) + bo,NPo,N{x), (3.1) 
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where a^^N > and bi^^jy G M are called the recurrence coefficients. By comparing 
(3.1) with the (unique) decomposition (1.18), one understands that the hypothesis 
(1.16) of Theorem 1.3 transposes to a sub-power growth condition as N ^ oo for 

max |afc,Ar| , max |^fc,Af| • (3-2) 

fceN: |-^-l|<£ ' fceN:|-|-l|<£ 

Here we shall combine our results with the work [41], where Kuijlaars and Van 
Assche obtain an explicit formula for the limiting zero distribution of general OPs for 
which the recurrence coefficients converge to some limit. More precisely, let us use 
the notation livcik/N^s (^k,N = c(s) when, for every sequences {kn)n£N and (A^n)neN 
such that kn,Nn — )• oo and kn/Nn — )• s as n — )■ oo, we have linin^oo Ck„,N„ = c(s). 
Let us also introduce the equilibrium measure of the interval [a,/3]. 



:,/3](d2;) 



' I lla£]ix)dx 

, for a < p, 

vr ^(i3-x){x-a) (3.3) 

6a if a = 13. 



Then the following holds. 



Theorem 3.1. Assume there exists e > and two continuous functions a{s), b{s) 
on [0, 1 + e) such that 

lim ak]\f = a{s), lim bkN = b{s), sG[0, 1 + e). (3.4) 

k/N^s ' k/N^s 

Then, W-almost surely, the empirical measure fi^ associated to the sequence of OP 
Ensembles converges weakly towards the probability measure 

The latter measure is defined in the obvious way, that is it evaluates any bounded 
and continuous function / on M to 



f{x)w\,, , , ^l (dx) ds. 

^ ' \b{s)-2a{s) ,b(s)+2a{s)\^ 'j 



As we shall observe below on several examples, all the OP Ensembles associated 
to classical OPs satisfy the conditions of Theorem 3.1, and one can recover the 
classical limiting distributions (semi-circle law, Marchenko-Pastur law, arcsine law 
...) from the formula (3.5). It is in that sense we claimed Theorem 3.1 provides a 
unified way to describe the almost sure convergence for classical OP Ensembles. 

Proof of Theorem 3.1. [41, Theorem 1.4] provides the weak convergence of z/^v to- 
wards (3.5). Let us show that i^n moreover converges to (3.5) in moments. Lideed, 
one can see from (3.1) that the spectral radius pN of 'kn^t^n-, which is the one of 
the matrix 

r 1 ^"-^ 

{xPk,N,Pm,N)L2(fj.,^) , , 

L J ft, 771=0 
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satisfies 

PN sup \ak,N\ + sup \hk,N\- 
k<N k<N 

Then, (3.4) yields sup^ /jat < +00, and Proposition 2.3 implies that the supports 
of the i^Ar's are uniformly bounded, from which the convergence in moments of i^at 
follows. Since (3.4) moreover provides that the sequences (3.2) are bounded, and 
because (3.5) has compact support, Theorem 3.1 follows from Corollary 1.5. □ 

We now provide a (non-exhaustive) list of recurrence coefficients for several 
rescaled classical OPs, from which one can check that Theorem 3.1 applies. The 
following formulas are easily obtained from the OP literature, see e.g. [34, Section 
9], and obvious change of variables. 





Orthogonal polynomial 


/iAf(dx) 


Parameters 


1 


Hermite 




none 


2 


Laguerre 


x^-e-^n[o,+oo)(x)dx 


a> -I 


3 


Jacobi 


(1 - a;)^"(l + a;)^^l[_i^i] {x)dx 


a,i3 > 


4 


Charlier 




a> 


5 


Meixner 




a £ (0,1), /3 > 





OP Ensemble 


(a/c.Af)^ 


bk,N 


1 


CUE 


k 
N 





2 


Wishart 




^+« 


3 


Random projections 






(2^+a+/3)2((2^+a+/3)2-^) 


(2^+a+/3)(2i±l+a+/3) 


4 


Longest increasing subsequence 




a + jf 


5 


Random Young diagrams 




T^(# + «(# + /3)) 



Note that these recurrence coefficients may still satisfy (3.4) if one lets the pa- 
rameters a, (3 depend on N in an appropriate way; e.g. letting a or /3 going to zero 
as N ^ 00. 

For more details concerning these OP Ensembles, we refer to [13] for the connec- 
tion between the product of random projections and Jacobi polynomials, to [33] for 
the problem of the longest increasing subsequence and Charlier polynomials, and to 
[32, 12] for the random Young diagrams and Meixner polynomials. 

It is finally easy to recover from Theorem 3.1 and the latter list of recurrence 
coefficients the almost sure convergence results for classical OP Ensembles. See also 
[41] for further computational examples. 
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4 Application to multiple orthogonal polynomials 



MOPs have been introduced in the context of the Hermite-Pade approximation of 
Stieltjes functions, which was itself first motivated by number theory after Hermite's 
proof of the transcendence of e, or Apery's proof of the irrationahty of ({2) and C(3), 
see [47] for a survey. For our purpose here, we will focus on the so-called type II 
MOPs, for which the zeros are of important interest since they are the poles of 
the rational approximants provided by the Hermite-Pade theory. These polynomi- 
als generalize orthogonal polynomials in the sense that we consider more than one 
measure of orthogonalization, and a class of classical MOPs such as multiple ver- 
sions of the Hermite, Laguerre, Jacobi, Charlier, Meixner, etc, polynomials emerged 
[14, 4, 3]. They are already the subject of many works where they are studied as 
special functions; we refer to the monograph [30] for further information. 

It turns out that even for the multiple Hermite or multiple Laguerre polynomials, 
no general description of the limiting zero distribution seems yet available in the 
literature. Our purpose in this section is to obtain such descriptions as a combination 
of our results with ingredients taken from free probability theory. 

Let us first introduce MOPs. 

4.1 Multiple orthogonal polynomials 

Let n he a Borel measure on M with infinite support and having all its moments. 
Consider r > 1 pairwise distinct functions wi, . . . ,Wr in Lp'{^). 

Definition 4.1. Given a multi-index n = (ni,...,nr-) G N^, the n-th (type II) 
MOP associated to the weights wi,. . . ,Wr and the measure ^ is the unique monic 
polynomial of degree ni + • • • + which satisfies the orthogonality relations 

J x^Pnix)wiix)fi{dx) = 0, 0<k<ni-l, 

■■ ■■ (4.1) 

j x^Pnix)wrix)n{dx) = 0, < A; < n,. - 1. 

Note that the existence/uniqueness of the n-th MOP is not automatic, and 
depends on whether the system of linear equations (4.1) admits a unique solution. 
We say that a multi-index n is normal if it is indeed the case. Since by taking r = 1 
we clearly recover OPs, we shall assume r > 2 in what follows. 

Let (n(^))ivgN = {n[^\ . . . ,ni^'')N^fi be a sequence of normal multi-indices 
which satisfies the following path-like structure. 

(a) For every G N, 

N 
i=l 

(b) For every € N and 1 < i < r, 

(Af+l) ^ {N) 



17 



(c) There exists i? G N such that for any G N and 1 < i < r, 

> ^(N) + 1. (4.2) 

(d) For every I < i < r, there exist qi, . . . ,qr G (0, 1) such that 

hm = q,. (4.3) 

We then write for convenience 

Piv(x) = P„(iv,(x), iVGN, (4.4) 

and observe that has degree A^. A question of interest is then to describe the 
weak convergence of the zero counting probabihty measure i^at of P/v as A'" — )• oo, 
defined as in (1.4) with zi, . . . , ztv the zeros of Pn{x), maybe up to a rescahng of 
the zeros. Before showing how our results answer that question in the case of the 
multiple Hermite and multiple Laguerre polynomials, we first need to introduce a 
few ingredients from free probability theory. 

4.2 Elements of free probability 

Free probability deals with non-commutative random variables which are indepen- 
dent in an algebraic sense. It has been introduced by Voiculescu for the purpose of 
solving operator algebra problems. We now just provide the few elements of free 
probability needed for the purpose of this work, and refer to [49, 1] for comprehensive 
introductions. 

For a probability measure A on R with compact support, let Kx be the inverse, 
for the composition of formal series, of the Cauchy-Stieltjes transform 

G.iz) = [ ^ (4.5) 



A:=0 ^ ^ 

and set the i?-transform of A by 

Rx{z) = Kx{z) - -. (4.6) 
z 

Definition 4.2. Let A and rj be two probability measures on M with compact sup- 
port. The free additive convolution of A and r/, denoted by A ffl ry, is the unique 
probability measure (on M with compact support) which satisfies 

RxBir,{z) = Rx{z) + R^iz). (4.7) 

Consider a probability measure A on [0, +oo) with compact support different 
from 6o. If x\ is the inverse for the composition of formal series of 



i^-A 0) - 1 = E ( / ^'^(d-)) A (4. 
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we then define the S'-transform of A by 

Sxiz) = ^xxiz). (4.9) 

Definition 4.3. Let A and r] be two probability measures on [0, +oo) with compact 
support and both different from 6o- The free multiplicative convolution of A and rj, 
denoted AKI??, is the unique probability measure (on [0, +oo) with compact support 
and different from 5q) which satisfies 

Sxm^iz) = Sx{z)Sr,{z). (4.10) 

For this work, the importance of the free additive and multiplicative convolutions 
relies on the following results due to Voiculescu, extracted from [1], which describe 
the limiting eigenvalue distribution of perturbed GUE and Wishart matrices. A 
random matrix X^v is distributed according to GUE(A^) if it is drawn from the 
space T-L]\f{C) N x N Hermitian matrices according to the probability distribution 



^ exp I - iVTr(X^)/2|dXAr, (4.11) 



Zn 

where dX^r stands for the Lebesgue measure on 7ii\f{C) ~ and Zjy is a nor- 
malization constant. It is said to be distributed according to WishartQ,(A^), where 
a > if a real parameter, if the probability distribution reads instead 

^det(X^)^-exp{ -iVTr(X;v)}l|x^>o}dXjv, (4.12) 

where Xjv > means that Xjv is positive semi-definite. The semi-circle distribution 
is defined by 

^sc(dx) = ^\/4-x2 l[_2,2](a;)dx, (4.13) 
and the Marchenko-Pastur distribution of parameter a > by 

^MP(a)(da;) = max(l — Q,0)(^o + min(Q;, 1)/Q,(x)dx, (4-14) 
where, with a± = {I i: -v/a)^) we introduced 

fa{x) = ^:^V («+ - - (4-15) 

Then the following holds. 

Theorem 4.4. Consider a sequence of uniformly bounded deterministic matrices 
(Ajv)jv, were Ajv is an N x N Hermitian matrix, and assume there exists a proba- 
bility measure A on M with compact support such that for all i gN, 



lim ^Tr (Ajv)^ = / x^A(dx) 

!>00 N J 



(a) If (X.]\[)n- is a sequence of independent random matrices with Xjv distributed 
according to GUE(A^), then for all ^ G N, 



lim — E 



Tr (Xat + An j 



x^fisc ffl A(dx). 
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(b) If (X.m)n is a sequence of independent random matrices with Xjv distributed 
according to WishartQ(A^), and if the Ajv 's are moreover positive semi-definite 
with A 7^ (5o, then for all i €zN, 



lim — E 

Af-5-oo N 



Tr (A 



N 



I- 



Vmph ^ A(dx) 



We are now in position to state the results of this section. 



4.3 Multiple Hermite polynomials 

Recall that if stands for the A^-th Hermite polynomial, that is the OP asso- 
ciated to /i(dx) = /^drc, then the zero counting probability distribution of 
its rescaled version Hn{^/Nx) is known to converge weakly towards the semi-circle 
distribution (4.13). 

Given r > 2 pairwise distinct real numbers ai, . . . , a,., consider the measure and 
the weights given by 

;u(dx) = e"^^/^dx, Wj{x) = e°-^'^, I < j < r. 

The associated MOPs are called multiple Hermite polynomials. For a sequence 
of multi- indices {n^^'>)N satisfying the path-like structure described in Section 4.1, 
denote by .■■■>'*'■) ^j^^ associated MOP as in (4.4). We shall prove the following. 

Theorem 4.5. Let vj^ he the zero probability distribution of the rescaled multiple 
Hermite polynomial 

i/^/^'^l•■■■•^"'■)(^/iVx). 
Then vj^ converges weakly as N ^ oo towards 

r 

^scffl 

J = l 

Although we introduced the i?-transform of a probability measure as a formal 
series, it is actually possible to define it as a proper analytic function, provided 
one restricts oneself to appropriate subdomains of the complex plane, and equal- 
ity (4.7) continues to hold, see [6, Section 5]. Then, since i?^gp(z) = z and the 
Cauchy-Stieltjes transform of Yll=i Qi^ai is explicit, one can obtain from (4.7) that 
the Cauchy-Stiejles transform G of fisc S ( Yl^j=i Qj^aj) is an algebraic function, by 
performing similar manipulations than in the proof of [7, Lemma 1] and concluding 
by analytic continuation. More precisely, one obtains that G satisfies the algebraic 
equation 

P{z,G{z))=0, zeC, (4.16) 
where the bivariate polynomial P{z, w) is given by 

r r r 

P{z,w) = wY[iz - w - ai) -'^Qi (z-w-ai). (4.17) 

1=1 i=l j=l,j^i 
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Probability measures for which the Cauchy-Stieltjes transform is algebraic have in- 
teresting regularity properties, see [2, Section 2.8], and are moreover suitable for 
numerical evaluation, see e.g. [25]. 

We now turn to multiple Laguerre polynomials, for which we provide a similar 
analysis. 



4.4 Multiple Laguerre polynomials 

If L^'* stands for the N-th Laguerre polynomial of parameter a > 0, that is the OP 
associated to fi{dx) = x"e~^l[o _|_oo)(x)dx, then it is known that the zero probability 

distribution i/jy of (Nx) converges weakly as — )• oo towards the Marchenko- 
Pastur distribution of parameter a (4.14)-(4.15). 

There exist two different definitions for the multiple Laguerre polynomials in the 
literature, see [30, Section 23.4]. We consider here the so-called multiple Laguerre 
polynomials of the second kind, which are defined as follows. Given r > 2 pairwise 
distinct positive numbers ai, . . . , and a > 0, consider 

Mdx) = x-e-n[o,+oo)(^)drE, wj{x) = e'-'+''^^^ , 1 < i < r, 
and, given a sequence of multi-indices (n(^))jv satisfying the path-like structure 



described previously, let ^(^■'^i '•••''^'■) |-,g i^j^g associated MOP as in (4.4). 



Theorem 4.6. Let vjsi he the zero probability distribution of the resettled multiple 
Lttguerre polynomittl 

j^(Na;Nai,...,Nar)^^^ 

Then vjsi converges wettkly as N ^ oo towards 

r 

i=i 

As it was the case for the ii-transform, the 5-transform can be defined as an 
analytic function, and (4.10) also holds on subdomains of the complex plane, see 
[6, Section 6]. Then, because 'S'^mp(c«) (-^) = (1 + ctz)~^ , one can also obtain from 
(4.10), taking care of the definition domains, that the Cauchy-Stieltjes transform G 
of /xmp(q) ^ ( J2j=i Qj^i/aj) satisfies the algebraic equation 

P{z,G{z))=0, zeC, (4.18) 

where P{z,w) is given by 
P{z,w) = wY\_{z (l-a + azw)) -'^Qi [z {1 - a + azw)) . (4.19) 

i=l * i=l j=l,j^i * 



4.5 Proofs 

Before providing proofs for Theorems 4.5 and 4.6, we first precise a few points 
concerning MOP Ensembles, that we introduced in Section 1.3. 
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A sequence of measures {^ly)^, weights wj^jy G L'^ifJ-N), 1 < i < '''5 and a 
path-hke sequence of multi-indices {n^^^)]\f induce a sequence of MOP Ensembles. 
Namely, for each one can associate random variables xi, . . . ,xn- distributed ac- 
cording to (1.7) where we chose for the multi-index n = n^^-*. For the monic 
polynomial Pk^N of degree k appearing in the kernel (1.9), one can choose the n^^^- 
th (type II) MOP associated with fijy and the wj^^^s. The associated biorthogonal 
functions Qk,N^s can then be constructed from the type I MOPs, see [30, Theorem 
23.1.6], and Assumption 1.1 is satisfied with q^v = 1- We moreover recall that the 
average characteristic polynomial xn equals Pn.n- 

In order to obtain growth estimates for the {xPk,NiQm,N) L^{fj,j^)^s, we now de- 
scribe a connection with the so-called nearest neighbors recurrence coefficients, which 
are in practice easier to compute. 



4.5.1 Nearest neighbors recurrence coefficients 

Van Assche [48] established for general MOPs, says associated to a measure /i and 
weights Wj's, that for every normal multi-index n there exist real numbers (a^f^)i<d<r 
and (6n^)i<d<r satisfying 



cPn(x) = Pn+e, + ak''Pn(x) + b'n^ P, 



(4.20) 



r 

vPnix) = Pn+e,, + a^n^ Pn{x) + ^ 6^^^ P„-e, (^ 



d=l 



where 



erf = ( 0, . . . , 0, 1, 0, . . . , 0) G W, l<d<r. 

d-l 

Note that this provides 

P„+e,(x)-P„+e^.(x) = (aL''^-aL'Vn(x), l<i,j<r. (4.21) 

With the path-like sequence of multi-indices (n'^'^))^^^ and allowing the u^j's and 
fi to depend on a parameter N, we write for convenience 

Then the following holds. 

Lemma 4.7. // there exists e > such that for every 1 < d < r the sequences 

{max max|a*-'jL .,( > , < max l^l'^Ll ^ , (4.22) 
are bounded, then so is the sequence 

{ , max \{xPk,N,Qm,N)L^(uN)\ } 

[k,mm:\j,-l\<e,\^-l\<e J 
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Proof. First, as a consequence of (4.2) and [30, (23.1.7)], we have 

(xPfc,iV,Qm,Af)L2(/,jv) = 0' m<k-R. 
Define the sequence {ik)keN taking its values in {1, ... , r} by 



n 



(k+l) _ „(fc) 



n^"' + e 



m G N. 



For a fixed k, which may be chosen as large as we want, (4.20) yields 



xPk,N{x) = Pk+l,N{x) + 4*';^Pfc,Ar(x) + ^ ^'^iV^n 



('')-ea,N 



d=l 



Then, since (4.21) provides for any 1 < d < r and m large enough 



^n(-)-ed,iv(a^) = Pm-l,N{x) + (a^--! 

we obtain inductively with (4.24) that 



(«m-l) 



) ~ea,N "n(™-i) -e^.A^^ ^n{— l) -e^.A^ 



(4.23) 



(4.24) 



r{x), 



xPk,N{x) = Pk+l,N{x) + a^^^^Pk,N{x) + ( ^ 4'!]v)^fc-l,A^( 



k-2 



k-1 



+ E I EC n ("11' 

m=k—R \ d=l l=m+l 



d=l 



d^i) 



)] Pm,N{x) + Rk,Nix), (4.25) 



where Rk,N is a polynomial of degree at most k — R — 1. By comparing (4.25) with 
the (unique) decomposition (1.18) and (4.23), we obtain explicit formulas for the 
{xPk^Ni Qm,iv)'s in terms of the nearest neighbor recurrence coefficients, from which 
Lemma 4.7 easily follows. □ 

4.5.2 Proof of Theorem 4.5 

Proof. Associate to the multi-indices {n^^^)]\[^f^ the (uniformly bounded) sequence 
(AAr)7veN of diagonal matrices 



diag( ai. 



, . . . , ai , ... , u/j., . . . , 



AN) 



AN) 



On the one hand, let (X.j\f)]sr be a sequence of independent random matrices, with 
Xjv distributed according to GUE(A^). If jj^ stands for the empirical measure 
associated to the eigenvalues of Yjy = Xat + A^r, then Theorem 4.4 (a) and (4.3) 
provide for any £ G N 



hm E 



lim 

N-^oo 



^E[Tr(Xjv + Aiv) 

ffl(E«A.)(da^)- (4-26) 
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On the other hand, observe from (4.11) that the random matrix Yjv is distributed 
on 7ij\f{C) according to 

^ exp I - iVTr(Y^ - 2A^Y^)/2|dY^, (4.27) 

where Z'j^ is a new normahzation constant. By performing a spectral decomposition 
in (4.27), integrating out the eigenvectors and using a confluent version of the Harish- 
Chandra-Itzykson-Zuber formula, Bleher and Kuijlaars [9] obtained that the random 
eigenvalues of Y^r form a MOP Ensemble, see (1.7), associated to the A^-dependent 
weights and measure 

fiNidx) = e-^^'/2^x, Wj^N{x) = e^'^^^, 1 < j < r, (4.28) 

and the multi-index n*^^-*. The average characteristic polynomial xn for that MOP 
Ensemble then equals the associated n^-^^-th MOP, which is seen from a change of 

variable to be H^/^"' ' - ' ^^^^ {^Nx) , up to a multiplicative constant. The weights 
in (4.28) form an AT system, from which it follows that any multi-index is normal, 
and that xn has real zeros, cf. [30, Chapter 23]. One moreover obtains from [48, 
Section 5.2] and a change of variables explicit formulas for the nearest neighbors 
recurrence coefficients associated to (4.28), 

Vat - ~ ' n - (77-1, . . . ,nr). 

Thus, Theorem 4.5 follows from (4.26), Lemma 4.7 and Corollary 1.6. □ 



4.5.3 Proof of Theorem 4.6 

Proof. The proof follows the same spirit as the proof of Theorem 4.5. Introduce the 
sequence of (uniformly bounded) diagonal matrices 

An = diag( 1/ai, . . . , 1/ai, . . . , l/a^, . . . , 1/a^ ), 



AN) 

^1 



AN) 



and let (Kn)n be a sequence of independent random matrices, where X^r is dis- 
tributed according to Wisharta(A^). With (i^ the empirical measure of the eigen- 

1 /2 1 /2 

values of Y^r = -^n ^a^-'^tv > Theorem 4.4 (b) and (4.3) then provide for all £ E N 



lim E 



J r.N 



x'fi'^dx) 



/r 

^(E'?iV.)(dx). (4.29) 
i=i 

Now, observe from (4.12) that Y^r is distributed on ^{^{C) according to 

^det(Y;v)^"exp { - 7VTr(A]^iYjv) ^jdY^, (4.30) 

where Z'j^ is a new normalization constant. Similarly than for the Hermite case, the 
eigenvalues of Yat form a MOP Ensemble associated to 



^Ar(dx) = X^°'e ^^dx, Wj^N{x) 



(4.31) 
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and the multi- index n*^^), see [10]. The average characteristic polynomial xn is 
then the n(^)-th MOP associated to (4.31), which is lSJ""'^"' ' ■ '^"''^ (x) up to 
a multiplicative constant. The weights in (4.31) form an AT system so that any 
multi-index is normal and XN has real zeros. If we denotes |n| = ni + ■ ■ ■ + for 
n G N'', then one obtains from [48, Section 5.4] that the nearest neighbors recurrence 
coefficients for (4.28) read 

id) _ nrf(|n| + jVg) (rf) _ |n| +iVa + l nj 

Theorem 4.6 finally follows from (4.29), Lemma 4.7 and Corollary 1.6. □ 

Remark 4.8. Having in mind the proofs of Theorems 4.5 and 4.6, it would be of 
interest to find out if there exists a matrix model for the multiple version of the 
Jacobi polynomials, the Jacobi-Piheiro polynomials, which are related in a limiting 
case to the rational approximations of C,{k) and polylogarithms [30, Section 23.3.2], 
and then if it would be possible to describe its limiting zero distribution thanks to 
free convolutions. 
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